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Richard  D.  Peacock  and  Richard  A.  Dipert 


Abstract 


A computer  implementation  of  a model  to  predict  temperatures  on 
wall  and  wall  protector  surfaces  exposed  to  the  heating  of  an 
appliance  such  as  a solid  fuel  heating  appliance  is  described.  A 
steady  state  heat  transfer  model  with  flexibility  to  describe  a 
generalized  method  of  protection  for  a combustible  wall  surface  is 
presented  along  with  a computer  program  implementing  the  model. 

Good  agreement  was  found  comparing  the  model  predictions  with  data 
previously  collected  during  full  scale  experiments  conducted  to 
evaluate  the  effectiveness  of  generic  methods  of  wall  protection 
in  reducing  temperatures  on  combustible  wall  surfaces. 

Extensive  references  of  research  related  to  solid  fuel  heating 
safety  are  included. 


Key  words:  Chimneys;  fire  models;  fire  safety;  fire  tests;  flues;  heat 

transfer;  heating  equipment;  literature  reviews;  radiant  energy;  stoves;  wood. 


1.  INTRODUCTION 

The  U.S.  Consumer  Product  Safety  Commission  and  the  U.S.  Department  of 
Energy,  as  part  of  a program  to  investigate  safety  risks  involved  with  the  use 
of  solid  fuel  burning  appliances,  have  sponsored  experimental  research  at  the 
Center  for  Fire  Research  (CFR)  at  the  National  Bureau  of  Standards  (NBS)  to 
identify  hazards  associated  with  solid  fuel  heating.  The  studies  were  con- 
ducted to  provide  information  to  improve  safety  practices  for  the  use  of  the 
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appliances,  and  to  provide  data  upon  which  to  base  improved  codes  and 
standards . 

During  the  first  years  of  the  program,  an  accident  survey,  literature 
review,  and  codes  and  standards  analysis  were  performed  to  establish  accident 
patterns , to  determine  the  types  of  risks  involved  with  the  use  of  wood 
burning  appliances,  and  to  ascertain  the  adequacy  of  existing  codes  and 
standards  in  addressing  these  risks  [1-3]. 1 Overwhelmingly,  conditions 
related  to  installation,  operation,  and  maintenance  were  responsible  for  the 
fire  incidents  studied.  Only  a small  percentage  of  the  fires  was  attributed 
to  product  design  or  product  defects.  Thus,  safe  installation  and  use  of  wood 
burning  appliances  is  a critical  requirement  for  preventing  fire  accidents 
involving  the  equipment.  Much  of  the  criteria  for  the  installation  and  use  of 
wood  burning  appliances  are  based  upon  data  developed  over  forty  years  ago  and 
do  not  provide  information  on  materials  of  construction,  or  appliances 
available  in  the  current  market. 

The  present  program  at  CFR  includes  research  on: 

• clearances  to  combustibles  from  appliances  and  chimney 
connectors  [4] ; 

• methods  of  protection  to  allow  reduced  clearances  to  walls  and 
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Number  in  brackets  refer  to  literature  references  listed  in 
section  8 at  the  end  of  this  report. 


ceiling  surfaces  exposed  to  radiant  heating  by  appliances  and 
chimney  connectors  [5]; 

• temperatures  developed  in  and  around  fireplaces  with  and  without 
fireplace  inserts  installed  [6,7]; 

• intensity  and  duration  of  chimney  fires  in  factory-built  and 
masonry  chimneys  [8] ; 

• temperatures  on  combustible  material  surrounding  chimney  con- 
nectors passing  through  walls  and  / or  connecting  to  chimneys  [9] 
and; 

• prediction  of  temperatures  on  surfaces  of  combustible  walls 
exposed  to  heating  from  a typical  radiant  heating  appliance. 

This  report,  one  of  a series  of  reports  providing  information  from  the 
experimental  program  on  wood  burning  safety  at  NBS , presents  the  results  of 
the  development  of  a computer  based  implementation  of  a model  to  predict 
temperatures  on  a wall  surface  exposed  to  heating  from  a radiant  heating 
appliance  such  as  a wood  stove.  The  resulting  computer  program  allows  the 
user  to  specify  thermal  protection  to  reduce  temperatures  on  the  wall  surface 
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2.  REVIEW  OF  PREVIOUS  WORK 


2.1  Fire  Incidents  Involving  Wood  Burning  Appliances 


Recent  statistics  on  fires  and  injuries  related  to  wood  burning 
appliances  are  alarming: 


Change 

From 

Previous  Dollar 


Year 

Fires 

Year 

Deaths 

Loss 

1978 

66,800 

250 

$134  million 

1979 

70,700 

+14% 

210 

$175  million 

1980 

112,000 

+58% 

350 

1981 

130,100 

+16% 

290 

$265  million 

1982 

139,800 

+7% 

250 

$257  million 

1983 

140,600 

+0 . 6 % 

280 

$296  million 

Source : 

U.S.  Consumer 

Product  Safety 

Commission  [10,  11] 

There  were  more  fires  in  solid  fuel  burning  equipment,  and  a larger 
percentage  increase  over  previous  years,  than  were  reported  for  any  other  kind 
of  heating  equipment  --  including  gas,  electric,  and  oil  burning 
appliances  [10-11]. 

This  marked  increase  is  attributed  to  the  growing  installation  and  use 


of  wood  burning  stoves  in  homes  throughout  the  United  States  and  the  fact  that 
most  homes  are  made  of  combustible  construction.  Clearly,  accidental  fires 
from  wood  burning  appliances  are  an  increasingly  important  problem. 
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2.2  Clearances  In  Existing  Codes  and  Standards 

Recommendations  for  minimum  acceptable  clearances  to  combustible 
materials  for  the  installation  of  chimney,  chimney  connectors,  and  appliances 
are  specified  in  the  various  model  building  codes  and  recommended  practices 
manuals.  Reference  [12]  is  typical  of  the  specifications  found  in  the  codes. 
For  simplicity,  a single,  hopefully  conservative  clearance  is  given  for  each 
type  of  appliance  installed  without  protection.  No  allowance  is  made  for  the 
size,  heat  output,  heat  transfer  characteristics  or  other  features  unique  to 
individual  models.  Similarly,  only  a few,  specific  methods  of  protection 
employed  to  allow  reduction  of  these  clearances  are  recommended. 

Typically,  0.91  m of  clearance  is  specified  between  radiant  heaters  and 
unprotected  combustible  construction.  For  residential  solid  fuel  chimneys, 
typically  51  mm  of  clearance  is  required.  Chimney  connectors  for  solid  fuel 
burning  residential  appliances  require  a clearance  of  at  least  0.46  m to 
combustible  materials.  However,  as  with  appliances,  these  clearances  may  be 
reduced  by  the  use  of  appropriate  protection  applied  either  to  the  appliance 
or  to  the  combustible  surface. 

The  experimental  basis  for  these  code  requirements  is  not,  in  many 
cases,  quite  so  clear.  Several  experimental  studies  have  been  carried  out  to 
determine  minimum  acceptable  clearances  to  combustible  materials.  Voigt  [13], 
in  a 1933  publication,  recommends  a minimum  clearance  of  0.30  m for  chimney 
connectors  0.23  m in  diameter.  A more  extensive  study,  performed  by 
Underwriters  Laboratories  in  1943  [14] , presents  minimum  safe  clearances  for 
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both  unprotected  surfaces  and  surfaces  protected  by  various  methods. 

Distances  at  which  a maximum  temperature  rise  of  50° C above  room  temperature 
is  reached  are  presented  as  a function  of  the  temperature  of  the  exposed  face 
of  a heat  producing  appliance.  The  relative  protection  afforded  by  various 
materials  used  as  heat  barriers  between  the  appliance  and  combustible  surfaces 
is  also  examined.  Lawson,  Fox,  and  Webster  [15]  and  Lawson  and  Simms  [16] 
have  studied  the  heating  of  wall  panels  and  wood  by  radiation.  With 
experimentation  and  theoretical  predictions,  they  present  safe  clearances 
between  flue  pipes  and  wall  surfaces  as  a function  of  the  pipe  diameter  and 
the  pipe  surface  temperature.  To  maintain  a maximum  wall  temperature  of 
100°C,  0.15  m pipe  should  not  exceed  350°C  in  surface  temperature  at  a 
clearance  of  0.46  m [15].  Loftus  and  Peacock  [5]  present  the  results  of 
research  studying  clearances  and  methods  of  protection  for  wall  and  ceiling 
surfaces  exposed  to  radiant  heating  appliances.  A number  of  recommended 
methods  of  protection  to  reduce  temperatures  on  combustible  wall  and  ceiling 
surfaces  to  acceptable  levels  were  found.  In  the  study,  appliance  surface 
temperatures  from  300  to  450°C  were  used. 

These  experimental  studies  established  limits  for  two  important  para- 
meters: appliance  surface  temperature  and  clearance  to  combustibles  for 
unprotected  and  protected  surfaces.  Maximum  appliance  surface  temperatures 
for  the  appliances  studied  ranged  from  300  to  450° C;  average  appliance  surface 
temperatures  from  200  to  250°C.  Minimum  safe  wall  clearances  for  unprotected 
surfaces  ranged  from  0.31  to  0.91  m.  Most  of  the  current  code  provisions  are 
only  adequate  for  maximum  appliance  surface  temperatures  up  to  300  to  350° C. 
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2.3  Temperatures  Developed  in  Heating  Systems 

Tests  made  with  prefabricated  porcelain- enameled  metal  chimneys  for 
solid  or  liquid  fuel  furnaces  [17,18]  established  a limiting  temperature  rise 
of  190°C  on  the  outer  surface  of  the  chimney  for  a flue  gas  temperature  of 
537°C.  With  this  limitation,  wood  framing  space  51  mm  or  more  away  from  the 
chimney  was  considered  safe.  Satisfactory  insulation  of  the  chimneys  to 
reduce  the  outer  surface  temperatures  to  acceptable  levels  was  obtained  with 
asbestos  paper  plies  totalling  about  45  mm  in  thickness.  In  the  same  study, 
some  asbestos  cement  pipe  coverings  were  also  found  to  reduce  heat 
transmission  to  the  extent  required  for  safety  of  nearby  combustibles. 

To  establish  performance  requirements  for  lightweight  prefabricated 
chimneys,  tests  were  conducted  with  lined  and  unlined  masonry  chimneys  having 
102  mm  thick  walls  [19,20].  Hazardous  conditions  on  wood  framing  spaced  51  mm 
away  from  the  chimney  were  noted  with  a continuous  flue  gas  temperature  of 
482° C for  the  unlined  chimney  and  592° C for  the  lined  chimney.  However,  these 
hazardous  conditions  were  not  reached  in  the  lined  chimney  tests  until  after 
13  hours.  In  order  to  study  operating  conditions  with  typical  fuels,  a number 
of  firing  tests  [21]  were  conducted  with  heating  appliances  known  to  give  high 
flue  gas  temperatures,  using  wood  and  soft  coals  as  fuels.  With  a coal- 
fired,  jacketed  type  heater,  gas  temperatures  ranging  from  648  to  704°C  were 
measured  for  an  hour  or  more  in  the  flue  at  the  ceiling  level  above  the 


heater . 
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Lawson,  Fox,  and  Webster  [15]  presented  results  of  tests  to  measure 
surface  temperature  of  flue  pipes.  Measured  for  a variety  of  flue  systems 
using  solid  fuels  --  mostly  coal  and  coke  --  they  report  temperatures  of  about 
150°C  under  "normal"  conditions  and  temperatures  as  high  as  815°C  for  over 
fire  conditions. 

Fox  and  Whittaker  [21]  report  temperatures  on  metal  flues  of  several 
heating  appliances  operated  over  a range  likely  to  encountered  in  normal  use. 
Maximum  flue  pipe  surface  temperatures  ranged  from  704  to  815 °C  at  the  appli- 
ance flue  outlet,  360  to  510° C at  a distance  of  0.91  m from  the  appliance  flue 
outlet,  and  287  to  326°C  at  a distance  of  1.8  m from  the  appliance  flue 
outlet . 

Shoub  [17]  concluded  that  combustible  materials  will  be  ignited  if 
maintained  in  continued  contact  with  a masonry  chimney  of  120  mm  wall 
thickness  with  flue  gas  temperatures  of  400° C. 

In  tests  for  the  Department  of  Energy  [4] , temperatures  ranging  from  297 
to  436°C  during  normal  operation  and  377  to  693°C  during  over  fire  conditions 
were  noted  on  the  surfaces  of  several  wood  burning  appliances  when  tested  to 
prescribed  test  methods  [22].  A total  of  11  different  short  term  tests, 
ranging  from  1.9  to  25.6  hours  duration,  were  conducted  to  establish  normal 
firing  conditions  in  wood  burning  appliances  [23].  An  examination  of  the 
data  from  these  tests  shows  spikes  occurring  at  the  beginning  and  end  of  tests 
and,  apparently,  whenever  the  door  to  the  stove  was  opened.  These  sharp  in- 
creases in  temperature  were  attributed  to  a "high  fire"  in  the  morning  and  to 
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the  rapid  increase  in  active  flaming  when  the  door  was  opened  for  refueling 
the  fire.  The  average  stove  surface  temperature  rise  for  normal  burning 
ranged  from  177  to  218°C,  flue  gases  from  140  to  269°C,  inner  chimney  wall 
surface  119  to  241° C,  and  the  outer  chimney  wall  surface  14  to  48° C. 


2.4  Limiting  Safe  Temperatures  on  Combustible  Surfaces 

Listings  of  heat  producing  appliances  and  methods  for  setting  clearances 
between  appliances  and  combustible  surfaces  are  based  upon  Underwriters 
Laboratories  listings  [22]: 

• maximum  temperature  rise  of  65 °C  above  room  temperature  on  exposed 
surfaces;  and 

• maximum  temperature  rise  of  50° C above  room  temperature  on  unexposed 
surfaces,  such  as  beneath  the  appliance,  floor  protector,  or  wall 
mounted  protective  device. 

These  requirements  are  based  upon  the  fact  that  while  the  ignition 
temperature  of  wood  products  is  generally  quoted  to  be  on  the  order  of  200° C 
[24],  wood  that  is  exposed  to  constant  heating  over  a period  of  time  may 
undergo  a chemical  change  resulting  in  a much  lowered  ignition  temperature  and 
increased  potential  for  self - ignition. 

Mitchell  [25]  presents  data  on  wood  fiberboard  exposed  to  temperatures 
as  low  as  109°C  that  resulted  in  ignition  after  prolonged  exposure.  MacLean 


10 


[26,27]  reports  charring  of  wood  samples  at  temperatures  as  low  as  93° C.  He 
concludes  that  wood  should  not  be  exposed  to  temperatures  appreciably  higher 
than  66 °C  for  long  periods.  McGuire  [28]  suggests  that  the  maximum  safe 
temperatures  on  the  surface  of  a combustible  material  adjacent  to  a constant 
heat  source  should  be  no  more  than  100° C. 

Clearly,  the  ignition  of  wood  at  moderately  elevated  temperatures  is  a 
complex  phenomenon;  the  time  of  exposure  is  indeed  an  important  parameter 
[29,30].  While  exact  limits  recommended  in  the  literature  vary  due  to 
exposure  time  and  details  of  the  tests  conducted,  the  numerous  documented 
fires  involving  the  ignition  of  wood  members  near  low  pressure  steam  pipes 
[31]  suggest  an  upper  temperature  limit  for  wood  exposed  to  long-term  low- 
level  heating  should  not  be  appreciably  higher  than  100°C. 

Nearby  combustible  materials  other  than  wall  and  ceiling  surfaces,  such 
as  chairs  or  draperies  must  also  be  kept  a sufficient  distance  from 
combustible  materials  to  prevent  ignition  of  the  materials.  The  testing 
standards  and  model  codes  treat  all  combustibles  with  the  same  requirements. 
Thus,  the  0.91  m clearance  requirements  in  NFPA  211  and  maximum  temperature 
rise  requirements  in  the  Underwriters  Laboratories  testing  standards  apply 
equally  well  to  other  combustibles  as  well.  Similarly,  a temperature  limit  of 
100°C  is  more  than  adequate  to  protect  most  combustibles  used  in  furnishings. 
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2.5  Data  for  Model  Validation 

Much  of  the  available  literature  related  to  wood  heating  safety  provides 
a significant  amount  of  data  that  can  be  used  to  compare  theoretical 
predictions  with  experimental  measurements.  In  the  above  literature, 
references  [4],  [5],  and  [9]  provide  measurements  of  appliance  surface 
temperatures  and  wall  surface  temperatures  and  provide  a full  description  of 
the  experimental  setup  to  be  modeled.  Reference  [22]  provides  acceptable 
limits  on  combustible  surface  temperatures  for  use  in  predicting  minimum 
clearances , maximum  appliance  surface  temperatures , minimum  acceptable 
protector  thermal  properties,  and  the  like.  These  references  will  provide  the 
majority  of  data  for  the  comparison  of  theoretical  predictions  with  measured 
temperatures  on  surfaces  of  appliances , on  protector  systems , and  on 
combustible  walls. 


3.  THEORETICAL  BASIS  FOR  THE  MODEL 

Figure  1 presents  a schematic  diagram  of  a heating  appliance  / wall 
system  with  an  arbitrary  protection  system  between  the  appliance  and  the 
wall.  Heat  transfers  from  the  hot  stove  surface  through  any  intervening 
protection  to  the  wall  surface,  through  the  wall,  and  to  the  cooler  surround- 
ings. A few  assumptions,  reasonable  to  the  system  being  modeled,  simplify  the 
model  considerably: 
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• The  stove  is  operating  at  steady  state  conditions  (thus,  we  assume 
the  stove  has  been  operating  for  a period  of  time  and  has  reached 
a steady  operating  condition) . 

• Stove  is  at  a constant  uniform  surface  temperature. 

• Heat  transfer  through  air  spaces  in  the  system  takes  place  by 
radiation  and  convection  only. 

• Heat  transfer  through  solids  in  the  system  takes  place  by  conduc- 
tion only. 

With  these  assumptions,  a one - dimens ional , steady  state  model  of  the 
stove  / protector  / wall  heat  transfer  is  appropriate.  The  only  loss  in 
generality  of  the  predictive  capability  of  the  model  is  the  inability  to 
predict  any  time  dependent  behavior  of  the  system.  Since  the  intended  purpose 
of  the  model  is  to  study  the  fire  safety  of  the  stove  / protector  / wall 
system  under  worst  case  conditions,  this  loss  is  acceptable.  By  assuming 
steady  state  conditions  with  a constant  stove  temperature,  the  worst  case 
conditions  will  be  modeled. 


3.1  Radiative  Heat  Transfer 


For  heat  exchange  between  two  surfaces,  the  net  radiative  heat  transfer 


between  surfaces  1 and  2 is  given  by  [32]: 
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F: 2 , the  configuration  factor  for  radiative  exchange  between  surface  1 and 
surface  2,  is  defined  as  the  fraction  of  the  radiation  leaving  surface  1 which 
is  intercepted  by  surface  2.  Compilations  of  configuration  factors  are 
available  in  the  literature  [33,34].  For  the  stove  / wall  protector  geometry, 
the  following  equations  are  appropriate.  The  configuration  factor  for  a 
differential  element  to  a plane  parallel  rectangle  with  the  normal  to  the 
element  passing  through  the  comer  of  the  rectangle  is  given  by  [33]: 


F = — 
dA  2n 


(1+X2)' 


( 1+Y2 ) 


tan 


-1 


(1+X2)' 


+ 


tan 


-1 


(1+Y2) 


2^ 


(2) 


where  X = (width  of  rectangle)  / (distance  from  rectangle  to  element)  and  Y = 
(height  of  rectangle)  / (distance  from  rectangle  to  element) . Since  the 
configuration  factor  for  a surface  equals  the  sum  of  configuration  factors  for 
any  subdivision  of  the  surface,  the  configuration  factor  for  any  point  (X^-.Y^ 
on  the  wall  (or  first  protector)  surface  can  be  defined  from  equation  (2)  as 
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For  radiant  heat  exchange  between  two  identical,  directly  opposed  rec- 
tangles (such  as  two  protective  surfaces  of  the  same  size),  the  configuration 
factor  is  given  by  [33] 
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(4) 


3.2  Convective  Heat  Transfer 


For  convective  heat  transfer  between  two  surfaces  separated  by  an  air 
space,  the  net  heat  exchange  by  convection  is  given  by 


q = h (Tx  - T2) 


(5) 


where  h is  the  convective  heat  transfer  coefficient.  For  free  convection  at 
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the  surface  of  a vertical  surface  (such  as  the  wall,  appliance,  or  protector 
surface),  h can  be  found  from  the  equations  [35] 


m h L 

NUL  " “IT 


0.387  Ra, 


1/6 


0.825  + 


E 


+ (0.492/Pr) 


9/16“] 


8/27 


(6) 


3.3  Conductive  Heat  Transfer 

For  solids,  the  net  heat  exchange  by  conduction  is  given  by 

i - -r  (Ti  - V (7> 

While  for  some  materials,  the  thermal  conductivity,  k,  is  a function  of 
temperature,  for  most  materials,  the  assumption  that  k is  a constant  leads  to 
inconsequential  loss  in  generality.  For  instance,  for  aluminum,  the  thermal 
conductivity  changes  by  only  25%  over  a temperature  range  from  -170°C  to  well 
over  2000°C 
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3.4  Solution  of  the  Equations  for  an  Arbitrary  Protection  System 

At  steady  state,  the  heat  transferred  from  the  appliance  to  the  outside, 
figure  2,  is  equal  to  the  heat  transferred  through  any  element  or  group  of 
elements  within  the  system,  or 


< Ts  - V 


1TOTAL 


^OTAL 


( Ti  - w 

R. 

l 


(8) 


Equations  (1),  (5),  and  (7)  can  be  expressed  as  in  equivalent  resistance  forms 


as 
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For  the  current  problem,  these  equations  can  be  combined  into  two  as 


R. 
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(Airspaces)  R.  = — a (T  2 + i ) (T  "+'t"7 

h + 1 - e 1 1 - e2 

— r-  + t-  + r- 

1 12  2 


(13) 


Thus,  the  solution  method  with  a given  appliance  temperature,  Ts , and  a 
given  temperature  of  the  surroundings,  T0 , begins  by  assuming  temperatures  for 
the  intermediate  surfaces,  calculating  individual  resistances  from  equations 
(12)  and  (13)  to  determine  the  total  resistance,  R^j^  , calculating  the  total 
heat  flow  rate  from  equation  (8) , and  comparing  the  total  heat  flow  rate  to 
the  individual  heats.  If  the  calculated  individual  rates  are  sufficiently 
close  to  the  calculated  total  heat  flow  rate,  the  steady  state  solution  has 
been  obtained.  If  any  of  the  individual  heats  are  different  from  the  total 
heat,  new  estimates  for  each  of  the  intermediate  temperatures  are  made  and  the 
process  is  repeated  until  sufficient  agreement  has  been  found.  Reference  [36] 
provides  a number  of  methods  for  making  best  guess  estimates  for  the  next 
iteration  toward  the  solution.  For  the  current  problem,  a simple  linear 
search  is  used  where  one  variable  at  a time  is  changed  until  a local  optimum 
is  found. 

A number  of  literature  sources  are  available  to  allow  determination  of  the 
physical  data  required  for  the  input  (for  instance,  for  the  thermal  conduc- 
tivity and  emissivity) . References  [32]  and  [33]  provide  extensive  listings 
of  thermal  and  material  properties  for  common  materials. 
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4.  COMPUTER  IMPLEMENTATION  OF  THE  MODEL 


Included  as  Appendix  A is  a listing  of  the  FORTRAN  (written  in  ANSI 
standard  FORTRAN  77)  program  implementing  the  model  as  described  above.  The 
general  structure  of  the  program  is  illustrated  in  figure  3. 


Program  input  (detailed  in  Appendix  B)  takes  the  form  of  different  key 
words  with  arguments  to  specify  values  which  depend  upon  the  key  word.  In 
most  cases,  the  order  of  the  key  words  is  unimportant.  A description  of  each 
of  the  input  key  words  and  values  which  go  on  the  same  line  are  presented 
below: 


STOVE 

AIRSPACE 

FOR 

XWALL 

YWALL 

PROTECTOR 


<height>  <width>  <emissivity>  <temperature> 

<thickness>  <emissivity> 

<variable>  — <lower>  <upper>  <increment> 

<x  position> 

<y  position> 

<thickness>  <height>  <width>  <k>  <emissivity>  <temperature> 


the  program  follows: 


A sample  input  for 


FOR: 

STOVE: 

AIRSPACE: 

ALUMINUM  SHEET: 
AIRSPACE: 

ALUMINUM  SHEET: 
AIRSPACE: 

GYPSUM  WALLBOARD: 
F/G  INSULATION: 
BRICK  OUTSIDE: 
AIRSPACE: 

END: 


TSTOVE  = 473.15  673.15  10. 
0.9  0.5  0.5  473.15 
0.91  0.9 

0.00254  3.0  3.0  177.  0.9 
0.0254  0.9 

0.00254  3.0  3.0  177.  0.9 
0.0254  0.9 

0.0127  3.0  3.0  0.17  0.9 
0.0916  3.0  3.0  0.038  0.9 
0.0916  3.0  3.0  0.72  0.9 
273.15 


19 


This  input  specifies  a series  of  calculations  to  be  done  for  a stove  tempera- 
ture ranging  from  200° C to  400° C ( 473.15  to  673.15  K in  the  input  above)  for 
a wall  protection  system  consisting  of  two  aluminum  sheets  spaced  25  mm 
(0.00254  m in  the  input  above)  apart  and  placed  25  mm  from  an  insulated 
outside  wall  of  a house.  A clearance  of  0.91  m from  the  appliance  to  the  first 
aluminum  sheet  is  specified.  A sample  output  from  the  program,  using  this 
sample  data  set,  is  presented  in  table  1.  Calculated  temperatures  on  all 
surfaces  from  the  stove  to  the  outdoors  are  shown  along  with  the  specified 
sizes  and  thermal  properties  of  the  materials  used  for  the  walls  and  protec- 
tors. Execution  time,  of  course,  depends  upon  the  computer  in  use.  On  a 
typical  desk  top  personal  computer,  execution  of  the  above  test  case  required 
less  than  1.5  minutes. 


5.  COMPARISON  OF  THE  MODEL  PREDICTIONS  WITH  EXPERIMENTAL  DATA 

Figure  4 presents  a comparison  of  calculated  temperatures  on  the 
surfaces  of  wall  protectors  and  on  the  surfaces  of  combustible  walls  with 
experimentally  measured  values  taken  from  references  [4]  and  [5].  To  assess 
the  predictive  capabilities  of  the  model,  a range  of  conditions  from  the 
experimental  studies  were  simulated.  Wall  materials  used  in  these 
experimental  studies  ranged  from  uninsulated  gypsum  wallboard  to  a fully 
insulated  stud  wall  with  a brick  facing  on  the  exterior.  A number  of 
different  wall  protection  methods  were  taken  from  reference  [5],  varying  from 
a simple  sheet  metal  protector  to  a sheet  metal  / insulation  board  / air  spacr 
composite  protector  or  a ventilated  brick  protector.  A simple  cross  plot  of 
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the  calculated  values  and  the  experimental  values  illustrates  the  agreement  of 
the  model's  predictions  with  experimental  data  obtained  in  earlier  studies. 
Agreement  of  the  calculated  values  with  the  experimentally  measured  values, 
stated  as 


( ( ^calculated  " Measured  ) / Calculated  ) * ^0 

with  temperatures  expressed  in  absolute,  averaged  within  less  than  1 percent. 
Individual  agreement,  however  ranged  from  5 percent  low  (calculated  values 
lower  than  experimental)  to  4 percent  high.  Much  of  the  disparity  in  the 
comparison  can  be  explained  by  the  choice  of  ambient  conditions  for  the 
experimental  tests.  All  the  data  in  the  two  reports  were  described  in  terms 
of  temperature  rise  above  ambient  conditions.  Since  the  experimental 
calculations  are  based  upon  absolute  temperatures,  some  assumptions  had  to  be 
made  for  the  ambient  temperatures  in  the  surroundings  during  the  tests.  A 
variation  in  ambient  temperature  of  15° C could  change  the  calculated  surface 
temperatures  on  the  wall  surfaces  by  as  much  as  ±10  percent.  Thus,  the 
agreement  illustrated  in  figure  5 is  excellent  in  light  of  the  possible 
variation  in  the  calculations  depending  upon  the  assumed  ambient  temperature. 

6.  MODEL  CAPABILITIES  AND  EXAMPLES 

A theoretical  model  for  predicting  the  heat  transfer  between  the 
appliance  and  the  wall  surfaces  can  be  a useful  tool  not  only  in  the  design  of 
appliances  and  wall  protection  devices  but  also  in  the  design  of  future 
experiments  to  study  clearances  and  reduced  clearances  for  wood  burning 


appliances.  This  section  presents  some  examples  of  the  use  of  the  model  in 
predicting  temperatures  and  clearances  from  combustibles  for  both  protected 
and  unprotected  wall  surfaces. 
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6.1  Heat  Transfer  From  Appliance  to  an  Unprotected  Wall 

Figure  5 shows  calculated  wall  surface  temperatures  as  a function  of 
appliance  / wall  clearance  for  a medium  size  appliance  (an  appliance  0.5  by 
0.5  m on  the  side  parallel  to  the  wall  surface)  adjacent  to  an  unprotected 
wall  surface  for  appliance  surface  temperatures  from  150  to  350° C.  For  these 
calculations,  the  outside  air  temperature  (temperature  of  the  surroundings) 
was  assumed  equal  to  0°C.  The  wall  consisted  of  12  mm  gypsum  wallboard,  a 92 
mm  stud  space  with  glass  fiber  insulation,  and  a 92  mm  common  brick  facing  on 
the  outside  of  the  wall  exposed  to  the  outdoors.  At  an  appliance  clearance  of 
0.91  m,  appliance  surface  temperatures  greater  than  about  300  °C  would  lead  to 
temperatures  on  the  wall  in  excess  of  the  recommended  limit  [22]  of  50° C above 
room  ambient  temperature  at  a point  on  the  wall  directly  centered  behind  the 
appliance.  Since,  in  an  earlier  study  [4],  average  appliance  surface 
temperatures  of  about  200  °C  were  noted,  a sufficient  margin  of  safety  is 
indicated  for  an  appliance  this  size. 
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6.2  Heat  Transfer  From  Appliance  to  a Sheet  Metal  Protected  Wall 

Figures  6 and  7 show  calculated  wall  surface  temperatures  as  a function 
of  appliance  / wall  clearance  for  a medium  size  appliance  (an  appliance  0.5  by 
0.5  m on  the  side  parallel  to  the  wall  surface)  adjacent  to  a protected  wall 
surface  for  appliance  surface  temperatures  from  150  to  350° C.  As  before,  the 
outside  air  temperature  was  assumed  equal  to  0°C.  The  wall  protector 
consisted  of  two  sheets  of  aluminum  (2.5  mm  in  thickness)  separated  by  a 
ventilated  25  mm  air  space.  The  wall  protector  was  spaced  from  the  wall  by  a 
ventilated  25  mm  air  space.  The  wall  consisted  of  12  mm  gypsum  wallboard,  a 
92  mm  stud  space  with  glass  fiber  insulation,  and  a 92  mm  common  brick  facing 
on  the  outside  of  the  wall  exposed  to  the  outdoors.  With  the  surfaces  of  the 
protector  painted  black  and  at  an  appliance  clearance  of  0.91  m,  appliance 
surface  temperatures  greater  than  about  300° C would  lead  to  temperatures  on 
the  wall  in  excess  of  the  recommended  limit  of  50° C above  room  ambient 
temperature  --  in  fact,  about  the  same  as  the  case  with  no  protection. 

However,  when  the  surfaces  of  the  protector  are  left  unpainted  (shiny  aluminum 
surfaces),  appliance  surface  temperatures  higher  than  350°C  are  required  to 
raise  the  temperature  of  the  wall  surface  above  acceptable  limits. 

Conversely,  the  clearance  of  the  appliance  to  the  wall  could  be  reduced  from 
0.91  m to  0.3  m with  an  average  appliance  surface  temperature  of  200° C. 
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6.3  Heat  Transfer  From  Appliance  to  a Masonry  Protected  Wall 

Figure  8 shows  calculated  wall  surface  temperatures  as  a function  of 
appliance  / wall  clearance  for  the  same  appliance  (an  appliance  0.5  by  0.5  m 
on  the  side  parallel  to  the  wall  surface)  adjacent  to  a protected  wall  surface 
for  appliance  surface  temperatures  from  150  to  350°C.  Again,  the  outside  air 
temperature  was  assumed  equal  to  0°C.  The  wall  protector  consisted  of  a 92  mm 
thick  solid  brick  wall  was  spaced  from  the  wall  by  a ventilated  25  mm  air 
space.  The  wall  consisted  of  12  mm  gypsum  wallboard,  a 92  mm  stud  space  with 
glass  fiber  insulation,  and  a 92  mm  common  brick  facing  on  the  outside  of  the 
wall  exposed  to  the  outdoors.  With  the  surfaces  of  the  protector  painted 
black  and  at  an  appliance  clearance  of  0.91  m,  appliance  surface  temperatures 
greater  than  about  300° C would  lead  to  temperatures  on  the  wall  in  excess  of 
the  recommended  limit  of  50°C  above  room  ambient  temperature  --  again  not 
significantly  lower  than  in  wall  surface  temperatures  than  for  the  unprotected 
wall.  Note,  however,  one  of  the  major  thermal  characteristics  of  a masonry 
wall  protection  system  --  high  thermal  mass  --  is  not  accounted  for  in  a 
steady  state  prediction. 


7.  USES  FOR  AND  LIMITATIONS  OF  THE  MODEL 

A model  was  developed  to  predict  temperatures  on  protected  and  unpro- 
tected wall  surfaces  exposed  to  heating  from  a (primarily)  radiant  heating 
appliance.  A one  - dimens ional , steady  state  model  of  appliance  / protector  / 
wall  heat  transfer  showed  agreement  within  an  average  of  less  than  1 perccmt 
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when  compared  to  experimental  results  from  earlier  laboratory  studies.  A 
range  of  building  materials  typical  of  residential  construction,  along  with  a 
number  of  different  wall  protection  methods  were  simulated  in  the  comparison. 
As  a guideline  to  the  range  of  applicability  of  the  model,  the  variations  of 
the  thermal  properties  used  in  the  comparisons  were 


k: 

0.038 

to 

177 

W/m 

K 

e : 

0.1 

to 

0.9 

airspace  thickness: 

0.1 

to 

1.0 

m 

solid  thickness: 

0.0025 

to 

0.1 

m 

Improvements  in  the  program  implementation  of  the  model  are  possible. 

The  simple  linear  search  for  the  solution  of  the  equations  is  certainly  not 
the  most  efficient  method.  A number  of  search  methods  have  been  described  in 
the  literature  [36].  An  n-dimensional  simplex  search  where  the  next  guess  for 
a given  variable  depends  upon  the  values  of  the  other  n-1  variables  would 
improve  the  execution  speed  of  the  program.  Additional  improvements  would  be 
realized  with  any  of  a number  of  acceleration  methods,  also  available  in  the 
literature  [36].  Of  course,  a more  complicated,  harder  to  understand  and 
modify  program  would  result. 
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Table  1.  Sample  Output  From  STOVE 

DOUBLE  ALUMINUM  PLATE,  ALL  SURFACES  PAINTED  BLACK 
NUMBER  OF  NODES  IN  CALCULATION:  10 

POINT  ON  WALL:  (X):  0.000  (Y) : 0.000 

SERIES  OF  CALCULATIONS  FOR  VARIABLE  TSTOVE:  473.150  673.150  20.000 


I 

MATERIAL 

HEIGHT 

(m) 

WIDTH 

(m) 

THICK 

(m) 

EMISS 

K 

(W/m-K) 

TEMPERATURE 

(°C) 

0 

STOVE 

0.50 

0.50 

0.90 

200.00 

1 

AIRSPACE 

0.91 

0.90 

2 

ALUMINUM  SHEET 

3.00 

3.00 

0.00 

0.90 

177.000 

3 

AIRSPACE 

0.03 

0.90 

4 

ALUMINUM  SHEET 

3.00 

3.00 

0.00 

0.90 

177.000 

5 

AIRSPACE 

0.03 

0.90 

6 

GYPSUM  WALLBOARD 

3.00 

3.00 

0.01 

0.90 

0.170 

7 

INSULATION 

3.00 

3.00 

0.09 

0.90 

0.038 

8 

BRICK  OUTSIDE 

3.00 

3.00 

0.09 

0.90 

0.720 

9 

AIRSPACE 

0.00 

1.00 

0.00 

0:  STOVE  / AIRSPACE 

1:  AIRSPACE  / ALUMINUM  SHEET 

2:  ALUMINUM  SHEET  / AIRSPACE 

3:  AIRSPACE  / ALUMINUM  SHEET 

4:  ALUMINUM  SHEET  / AIRSPACE 

5:  AIRSPACE  / GYPSUM  WALLBOARD 

6:  GYPSUM  WALLBOARD  / INSULATION 

7:  INSULATION  / BRICK  OUTSIDE 

8:  BRICK  OUTSIDE  / AIRSPACE 

9:  AIRSPACE 


TSTOVE 

( ' 

0) 

( 

1) 

( 

2) 

( 

3) 

( 

4) 

( 

5) 

( 

6) 

( 

7) 

( 8) 

( 9) 

(°, 

O 

(°< 

C) 

C' 

C) 

(° 

C) 

(° 

C) 

(°' 

C) 

C' 

C) 

(°< 

C) 

(°C) 

(°C) 

473. 

150 

200 

.0 

33 

.1 

33 

.1 

31 

.6 

31 

.6 

30 

.1 

29 

.3 

3 

.1 

1.7 

0.0 

493. 

150 

220 

.0 

38 

.7 

38 

.7 

37 

.0 

37 

.0 

35 

.3 

34 

.4 

3 

.6 

2.0 

0.0 

513. 

150 

240 

.0 

44 

.7 

44 

.7 

42 

.8 

42 

.8 

41 

.0 

39 

.9 

4 

.2 

2.3 

0.0 

533. 

,150 

260 

.0 

51 

.1 

51 

.1 

49 

.0 

49 

.0 

47 

.0 

45 

.7 

4 

.7 

2.6 

0.0 

553. 

,150 

280 

.0 

57 

.8 

57 

.8 

55 

.6 

55 

.6 

53 

.4 

51 

.9 

5 

.4 

2.9 

0.0 

573. 

,150 

300 

.0 

64 

.8 

64 

.8 

62 

.5 

62 

.5 

60 

.1 

58 

.5 

6 

.0 

3.2 

0.0 

593. 

,150 

320 

.0 

72 

.2 

72 

.2 

69 

.7 

69 

.7 

67 

.2 

65 

.4 

6 

.7 

3.6 

0.0 

613. 

,150 

340 

.0 

79 

.9 

79 

.9 

77 

.3 

77 

.3 

74 

.6 

72 

.6 

7 

.4 

3.9 

0.0 

633. 

,150 

360 

.0 

87 

.8 

87 

.8 

85 

.1 

85 

.1 

82 

.3 

80 

.1 

8 

.1 

4.3 

0.0 

653. 

,150 

380 

.0 

96 

.0 

96 

.0 

93 

.2 

93 

.2 

90 

.3 

87 

.8 

8 

.8 

4.7 

0.0 

673. 

,150 

400 

.0 

104 

.3 

104 

.3 

101 

.5 

101 

.5 

98 

.5 

95 

.8 

9 

.6 

5.0 

0.0 

Figure  1.  Heat  Transfer  From  Hot  Appliance  to  Cooler  Wall  Surface 
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igure  2.  Node  Model  Representation  of  Stove  / Protector  / Wall  System 
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Figure  4.  Comparison  of  Calculated  Wall  and  Protector  Surface  Temperatures  With  Experimental  Results 
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Figure  6 . Predicted  Wall  Surface  Temperatures  for  Heat  Transfer  From  Appliance  Surface  to  a Combustible 
Wall  Protected  With  a Double  Aluminum  Sheet  Wall  Protector  (Painted  Black) 
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Figure  7.  Predicted  Wall  Surface  Temperatures  for  Heat  Transfer  From  Appliance  Surface  to  a Combustible 
Wall  Protected  With  a Double  Aluminum  Sheet  Wall  Protector 
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Figure  8 . Predicted  Wall  Surface  Temperatures  for  Heat  Transfer  From  Appliance  Surface  to  a Combustible 
Wall  Protected  With  a Solid  Masonry  Wall  Protector  With  A Ventilated  Air  Space 
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Appendix  A: 


Program  Listing  of  STOVE 


1 

2 

3 

4 

5 

6 

7 

8 

g 

10 

li 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 


PROGRAM  STOVE 


C 


c 

@@@@@@@@ 

c 

@@@@  @ 

@ 

c 

@@@@ 

@@@@@@ 

@@@@@@ 

@@§  @ 

@@@@@@ 

c 

@@@@@@@@ 

@@@@ 

@@@  @ 

@@@  @ 

@ @@@ 

c 

@@@(s 

@@@@ 

@@@  @ 

@@@  @ 

@@@@@@ 

c 

@ @@@@ 

@@@@ 

@@@  @ 

@§@@ 

@ 

c 

@@@@@@@@ 

@@@@ 

@@@@@@ 

@@ 

@@@@@@ 

C 

C 

C THIS  PROGRAMS  CALCULATES  AN  ENERGY  BALANCE  ON  A STOVE  / WALL  SYSTEM 
C ASSUMING  STEADY  STATE  CONDITIONS,  WITH  AN  ARBITRARY  NUMBER  OF  WALL 
C PROTECTORS  BETWEEN  THE  STOVE  AND  THE  WALL 
C 

INCLUDE  ’ STOVE. CMN’ 

CHARACTER  CMET*5, MATL*60 ,VERSN*7 
LOGICAL  THRU1 , NOT YET , THEEND 
DATA  VERSN  /’86.0108’/ 

C 

C INITIALIZE  EVERYTHING  TO  DEFAULT  CONDITIONS,  READ  INPUT  FOR  CASE 
C 

QEPS=0. 00001 
EPS=0 . 0001 
TSTEP1=2. 

TSTEP2=2 . 

10  CALL  INPUT (THEEND) 

IF  (.NOT.  THEEND)  THEN 
XSTOVE=WSTOVE/2 . 

YSTOVE=0 . 

IF  ( . NOT . FULPRT ) THEN 

WRITE  (6,11)  (I ,FMATL(MATL(I) ,MATL(I+1) ) , 1=0, N-l) ,N,MATL(N) 
WRITE  (6,12)  VNAME, (1,1=0, N) 

WRITE  (6,8)  • ’,(’  CC)  ',1-0, N) 

WRITE  (6,8)  (’ ’ ,1— 1,N) 

END  IF 

NCALCS-MAX ( INT ( ( VUPPER-VLOWER ) / VINCR+0 . 5 ) +1 , 1 ) 

DO  100  ICALC=1 , NCALCS 
VARVAL=VLCWER+VINCR* ( ICALC- 1 ) 

C 

C SET  NEXT  VALUE  OF  VARIABLE  TO  BE  INCREMENTED 
C 

IF  (IVAR.EQ. 1)  THEN 
XWALL-VARVAL 
ELSE  IF  (IVAR.EQ. 2)  THEN 
YWALL=VARVAL 
ELSE  IF  (IVAR.EQ. 3)  THEN 
T ( 0 )=VARVAL 

ELSE  IF  (IVAR.EQ. 4)  THEN 
L ( I SUB ) -VARVAL 
IF  (1SUB.EQ. 1)  ZSTOVE-VARVAL 
ELSE  IF  (IVAR.EQ. 5)  THEN 
WIDTH ( ISUB ) =VARVAL 
IF  (ISUB.EQ.0)  WSTOVE=VARVAL 
XSTOVE=WSTOVE/2 . 

ELSE  IF  (IVAR.EQ. 6)  THEN 
HEIGHT ( ISUB )=VARVAL 
IF  (ISUB.EQ.0)  HSTOVE-VARVAL 
ELSE  IF  (IVAR.EQ. 7)  THEN 
K ( I SUB ) -VARVAL 
ELSE  IF  (IVAR.EQ. 8)  THEN 
ENNIS ( ISUB )=VARVAL 
END  IF 
C 

C INITIALIZE  THE  COUNTERS  TO  ZERO 
C 

ITER=0 

C 

C INITIALIZE  TEMPERATURE  GUESSES  TO  WORST  CASE  CONDITIONS 
C 


71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 


71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 


DO  20  1=1, N-l 
20  T(I )=T(N) 

TDIFF=T(0)-T(N) 

TGOAL=T(N) 

C 

C DETERMINE  NEXT  GUESS  FOR  TEMPERATURES 
C 

DO  30  1=1, N-l 
30  T( I )=T( I )+TDIFF/N 

C 

C CALCULATE  TOTAL  RESISTANCE  AND  HEAT  FLOW  FROM  GUESSED  TEMPERATURES 

C 

CALL  HEATS 
40  INNER=0 

C 

C CALCULATE  WITH  NEWLY  GUESSED  TEMPERATURES  AND  HEATS 
C 

DO  60  1=1, N 

TINC=ABS ( TDIFF / N/ TSTEP1 ) 

TINCI=TINC 

QSIGN=DSIGN ( 1 . DO , Q( I ) -QTOTAL ) 

C 

C GO  THROUGH  INNER  ITERATION  LOOP  AT  LEAST  ONCE  FOR  SOME  FORCED 
C IMPROVEMENT 
C 

THRU 1=. FALSE. 

50  IF  (.NOT.  (ABS( (Q( I )-QTOTAL) /QTOTAL) . LE. QEPS. AND. THRU 1) ) THEN 

INNER=INNER+1 
THRU1=.TRUE. 

IF  (ABS ( (Q( I ) -QTOTAL ) /QTOTAL ) . GT . QEPS ) NOTYET= . TRUE . 

IF  (MATL ( I ).NE. ’AIRSPACE’ ) THEN 
C 

C HEAT  TRANSFER  IS  CONDUCTION,  CALCULATE  TEMPERATURE 
C 

T ( I )=T ( I- 1 ) -QTOTAL*R( I ) 

ELSE 

C 

C HEAT  TRANSFER  IS  BY  CONVECTION  & RADIATION,  SEARCH  FOR  TEMPERATURE 
C 

T(I )=T(I)+QSIGN*TINC 
C 

C RECALCULATE  THE  HEAT  FLOW  AND  COMPARE  TO  TOTAL  HEAT 
C 

CALL  HEAT(I) 

IF  (DSIGN(l.D0,Q(I)-QTOTAL) .NE.QSIGN)  THEN 
TINCI=TINCI/TSTEP2 
TINC=TINCI 

QSIGN=DSIGN ( 1 . DO , Q( I ) -QTOTAL ) 

ELSE 

TINC=TINC*2. 

END  IF 
GO  TO  50 
END  IF 
END  IF 
C 

C RECALCULATE  TOTAL  HEAT  WITH  THE  NEW  TEMPERATURES  AND  PROCEED 
C 

CALL  HEATS 
60  CONTINUE 

C 

C IF  DEBUG  IS  ON,  PRINT  OUT  A SIM1ARY  OF  THE  ITERATION 
C 

ITER=ITER+1 

IF  (DEBUG. AND. NOT YET)  THEN 
WRITE  (6,7)  ITER, INNER 

WRITE  (6, ’ ( IX, A,F7 . 2) ' ) ’STOVE  TEMPERATURE:  ’ ,T(0)-273. 15 
WRITE  (6, ’ (1X,A,F8.2,A,F8.3)’ ) 'TOTAL  HEAT:  ’.QTOTAL, 

2 ’ TOTAL  RESISTANCE:  ’ , RTOTAL 

WRITE  (6,3) 

DO  70  1=1, N 


141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 


141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 


70 


C 

C 

C 

C 


80 


C 

C 

C 


90 


100 


1 

2 

3 

4 

5 

6 

7 

8 
9 
11 
12 


IF  (MATL(I).EQ. ’AIRSPACE’. AND. I. EQ.l)  THEN 
WRITE  (6,5)  MATL(I) ,T(I)-273.15,R(I),Q(I),CF(I) 

ELSE  IF  (MATL(I).EQ. ’AIRSPACE' .AND. I. NE.l)  THEN 

WRITE  (6,1)  MATL( I ) ,T(I )-273 . 15,R(I),Q(I),H(I),CF(I) 

ELSE 

WRITE  (6,2)  MATL(I) ,T(I)-273.15,R(I),Q(I) 

END  IF 

CONTINUE 

TDIFF=T(N)-TGOAL 

WRITE  (6, ’ (1X,A7,F10.3,A)’ ) ’ TDIFF:  ’ ,TDIFF, 

2 CMET(TDIFF/TGOAL,EPS) 

WRITE  (6,6)  'QMET:  ’ , (Q(I)-QTOTAL, 

2 CMET ( (Q ( I ) -QTOTAL ) / QTOTAL , EPS ) , 1=1 , N ) 

END  IF 

TDIFF=T(N)-TGOAL 

T(N)=TGOAL 

IF  ANY  OF  THE  INDIVIDUAL  HEATS  HAVEN’T  CONVERGED  ON  THE  FINAL  VALUE, 
CONTINUE  WITH  THE  PROCESS  UNTIL  ALL  ARE  EQUAL 

NOTYET=. FALSE. 

DO  80  1=1, N 

IF  (ABS((Q( I) -QTOTAL) /QTOTAL) .GT.EPS)  NOTYET= . TRUE . 

CONTINUE 

IF  (ABS (TDIFF )/TGOAL. GT.EPS)  NOTYET= . TRUE . 

IF  (NOTYET)  GO  TO  40 

THE  SOLUTION  HAS  BEEN  FOUND,  PRINT  IT  OUT 

IF  ( . NOT . FULPRT ) THEN 
WRITE  (6,9)  VARVAL, (T(I )-273 . 15, 1=0 ,N) 

ELSE 

WRITE  (6.4)  TITLE 

WRITE  (6, ’(IX, A, 13)’)  ’ITERATIONS:  ’.ITER 

WRITE  (6, ’ ( IX , 3A,F8 . 3 ) ’ ) 'VARIED  VARIABLE:  ’ , VNAME, ’ = ’, 

2 VARVAL 

WRITE  (6, ’ ( IX , A,F7 . 2) ' ) 'STOVE  TEMPERATURE:  ’ ,T(0)-273. 15 
WRITE  (6, ’ (1X,A,F8.2,A,F8.3)’ ) ’TOTAL  HEAT:  ’.QTOTAL, 

2 ’ TOTAL  RESISTANCE:  ’ , RTOTAL 

WRITE  (6,3) 

DO  90  1=1, N 

IF  (MATL(I).EQ. ’AIRSPACE’. AND. I. EQ.l)  THEN 

WRITE  (6,5)  MATL(I) ,T(I)-273.15,R(I),Q(I),Q(I) 

ELSE  IF  (MATL(I).EQ. ’AIRSPACE’ .AND. I. NE.l)  THEN 
QCONV=(T(I-l)-T(I))/CONV(I) 

QRAD=(T(I-1)-T(I))/RAD(I) 

WRITE  (6,1)  MATL(I),T(I)-273.15,R(I),Q(I),QCONV,QRAD 
ELSE 

WRITE  (6,2)  MATL(I) ,T(I)-273.15,R(I),Q(I) 

END  IF 
CONTINUE 

WRITE  (6,4)  'TEMPERATURES  ARE  OF  COOLER  SIDE  OF  NODE’ 

END  IF 
CONTINUE 
GO  TO  10 
END  IF 
STOP 

FORMAT  ( IX , A10 , 2X ,F10 . 2 , 3X, F10 . 3 , 2X.F10 . 3 , 3X, F10 . 3 , 3X ,F10 . 3 ) 

FORMAT  ( IX, A10 , 2X.F10 . 2 , 3X, F10 . 3 , 2X ,F10 . 3 ) 

FORMAT  ( ’ 0MATERIAL  TEMPERATURE  RESISTANCE  HEAT  CONV 

2ECTION  RADIATION ’ , / ) 

FORMAT  ( ’ 0’ ,A) 

FORMAT  ( IX, A10 , 2X ,F10 . 2 , 3X, F10 . 3 , 2X ,F10 . 3 , 16X ,F10 . 3 ) 

FORMAT  ( IX, A7 , 4 (F10 . 3 , A5, : ) , / , 10( 8X , 4 (F10 . 3 , A5 , :),/)) 

FORMAT  ( ’ OITERATION : ’,13,'  INNER  LOOP:  ’,13) 

FORMAT  (’  ’ ,A10,20(A7, : )) 

FORMAT  (’  ’ ,F10.3,20(F6.1,1X, :)) 

FORMAT  ( ’ 0’ , / , 20( ’ ’,16,’:  ’,A,/)/) 

FORMAT  (’  ’ , A10 , 20 ( ’ (’ ,12, ’)’,:) ) 

END 
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C 

C PURPOSE:  INPUTS  DATA  DESCRIBING  THE  STOVE  / WALL  PROTECTION  SYSTEM 

C TO  BE  MODELED 

C 

INCLUDE  ’ STOVE. CMN’ 

PARAMETER  (NVAR=8) 

CHARACTER  VAR(NVAR)*6 , IN*80 ,KEYWD*80 
INTEGER  VLEN(NVAR) 

LOGICAL  THEEND , HAVEX , HAVEY 

DATA  (VAR(I ) , VLEN(I ) , 1=1, NVAR)  /’XWALL  ’,5,’YWALL  ’ , 5 , ’ TSTOVE’ , 6, 
2 'THICK  ',5, 'WIDTH  ' , 5, 'HEIGHT' , 6 , 'K  ’,1,'EMMIS  ’,5/ 

C 

C LOOK  FOR  A KEYWORD  OR  END  OF  FILE 
C 

THEEND= . TRUE . 

HAVEX=. FALSE. 

HAVEY=. FALSE. 

TITLE=’ STOVE  / WALL  PROTECTOR  HEAT  TRANSFER  MODEL' 

N=0 
IVAR=0 
VNAME=’  ’ 

10  READ  (5, ’ (A80)’ ,END=30)  IN 
THEEND=. FALSE. 

IS=INDEX( IN, ' : ’) 

IF  (IS.GT.l)  THEN 
KEYWD=IN( 1: IS-1) 

IF  (KEYWD.NE. 'END' ) THEN 
C 

C XWALL  ...  X POSITION  OF  POINT  ON  THE  WALL 
C 

IF  (KEYWD.EQ. 'XWALL')  THEN 
ICHR=IS+2 

XWALL=FVAL ( IN , ICHR) 

HAVEX= . TRUE . 

C 

C YWALL  . . . Y POSITION  OF  POINT  ON  THE  WALL 
C 

ELSE  IF  (KEYWD.EQ. 'YWALL' ) THEN 
ICHR=IS+2 

YWALL=FVAL ( IN , ICHR) 

HAVEY=.TRUE. 

C 

C AIRSPACE  ...  AN  AIR  SPACE  BETWEEN  TWO  SOLID  PROTECTORS 
C 

ELSE  IF  (KEYWD.EQ. 'AIRSPACE' ) THEN 
N=N+1 

MATL(N)=' AIRSPACE’ 

ICHRt=IS+2 
V1=FVAL(IN, ICHR) 

V2=FVAL(IN, ICHR) 

V3=FVAL ( IN , ICHR ) 

IF  (V2.EQ.0. .AND.V3.EQ.0. ) THEN 
T (N)=V1 
EM1IS (N)=l. 0 
HEIGHT(N)=0. 

ELSE 
L (N)=V1 
EM1IS(N)=V2 

IF  (EFMIS(N).EQ.O.)  EPMIS(N)=1.0 
T (N)=V3 
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HEIGHT (N)=HEIGHT(N-1) 

WIDTH ( N ) =WIDTH ( N- 1 ) 

END  IF 
C 

C STOVE  . . . THE  HOT  STOVE  SURFACE 
C 

ELSE  IF  (KEYWD.EQ. ’STOVE’ ) THEN 
MATL(0)=’ STOVE’ 

ICHR=IS+2 

EMMIS  ( 0 ) =FVAL  ( IN , ICHR ) 

IF  (EFMIS(O).EQ.O.)  EMIS(N)=1. 0 
HEIGHT ( 0 )=FVAL (IN , ICHR) 

WIDTH ( 0 )=FVAL ( IN , ICHR ) 

T (0 )=FVAL(IN, ICHR) 

C 

C FOR  . . . SPECIFIES  A SERIES  OF  CALCULATIONS  TO  BE  DONE 
C 

ELSE  IF  (KEYWD.EQ. ’FOR’ ) THEN 
IVAR=0 
ISUB=0 

DO  20  1=1 , NVAR 

J=I NDEX ( IN , VAR ( I ) ( 1 : VLEN ( I ) ) ) 

IF  ( J .NE . 0 . AND . IN( J~1 : J-l) . EQ. ’ ’)  IVAR=I 
20  CONTINUE 

IR=0 

ICHR=IS+2 

IF  (IVAR.EQ. 0 ) THEN 
WRITE  (6,1)  IN 
GO  TO  10 
C 

C IF  IT'S  A VARIABLE  WITH  SUBSCRIPT,  MAKE  SURE  ONE’S  THERE 
C 

ELSE  IF  (IVAR.GE.4)  THEN 
IL=INDEX( IN , ’ ( ’ ) 

IR=INDEX( IN , ’ )’ ) 

IF  ( IL . EQ . 0 . OR . IR . EQ . 0 . OR . IL . GT . IR)  THEN 
WRITE  (6,2)  IN 
GO  TO  10 
END  IF 

ISUB=FVAL(IN, ICHR) 

END  IF 
ICHR=IR+1 

VLOWER=FVAL ( I N , ICHR ) 

VUPPER=FVAL ( IN , ICHR ) 

VINCR=FVAL( IN , ICHR) 

IF  ( (VLOWER . NE . VUPPER. AND . VINCR.EQ. 0 ) 

2 .OR. (VLOWER. EQ. VUPPER. AND. VINCR.NE.0) 

3 .OR. (VUPPER.GT.VLOWER.AND.VINCR.LT. 0) 

4 .OR. ( VUPPER.LT. VLOWER. AND. VINCR.GT.O))  THEN 

WRITE  (6,4)  IN 

GO  TO  10 
END  IF 

VNAME=VAR ( I VAR ) 

VNAME ( VLEN ( IVAR )+ 1 : VLEN ( I VAR ) +1+IR- IL )=IN ( IL : IR ) 
C 

C DEBUG  . . . TURN  DEBUG  PRINT  ON  OR  OFF 
C 

ELSE  IF  (KEYWD.EQ. ’DEBUG’ ) THEN 

IF  (INDEX (IN, ’OFF’ ) .NE. 0)  DEBUG=. FALSE. 

IF  (INDEX (IN, ’ON’ ) .NE.O)  DEBUG= . TRUE . 

C 

C PRINTOUT  . . . SPECIFY  LEVEL  OF  PRINTOUT 
C 

ELSE  IF  (KEYWD.EQ. ’PRINTOUT’ ) THEN 

IF  (INDEX(IN, ’FULL’ ) .NE.O)  FULPRT= . TRUE . 

IF  (INDEX(IN, ’FULL’ ).EQ.O)  FULPRT= . FALSE . 

C 

C TITLE  . . . SPECIFY  A TITLE  FOR  THE  PRINTOUT 
C 

ELSE  IF  (KEYWD.EQ. ’TITLE’ ) THEN 
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TITLE=IN ( IS+2 : 80 ) 

C 

C IF  IT’S  NOT  ONE  OF  THE  RECOGNIZED  KEYWORDS,  ASSUME  A SOLID  PROTECTOR 
C 

ELSE 
N=N+1 
ICHR=IS+2 

MATL(N)=IN( 1 : IS-1 ) 

L(N)=FVAL(IN, ICHR) 

HEIGHT(N)=FVAL(IN, ICHR) 

WIDTH ( N ) =FVAL ( I N , ICHR ) 

K(N)=FVAL(IN, ICHR) 

EffllS  ( N )=FVAL  ( IN , ICHR ) 

T (N)=FVAL(IN, ICHR) 

IF  (L(N).EQ.O. .OR.HEIGHT(N) .EQ.0. .OR. WIDTH(N) .EQ. 0 . .OR. 

2 K(N) . EQ. 0 . . OR.EMIIS(N) . EQ. 0 . ) THEN 

WRITE  (6,3)  IN 
N=N-1 
GO  TO  10 
END  IF 
END  IF 
GO  TO  10 
END  IF 
END  IF 
C 

C DATA  HAS  BEEN  READ  IN,  CHECK  CONSISTENCY  OF  DATA  INPUT 
C 

30  IF  (THEEND)  THEN 
RETURN 

ELSE  IF  (N.LE.l)  THEN 

WRITE  (6,*)  'DATA  INPUT  ERROR,  TOO  FEW  NODES  SPECIFIED.' 

ELSE 

IF  (DEBUG. AND. I VAR. NE.0)  THEN 
IF  (IVAR.LT.4)  THEN 

WRITE  (6,5)  VAR ( I VAR ) , VLOWER , VUPPER , VINCR 
ELSE 

WRITE  (6,6)  VAR(IVAR),ISUB, VLOWER, VUPPER, VINCR 
END  IF 
END  IF 
DO  40  1=1, N 

IF  (DEBUG)  CALL  PSUM(I) 

IF  (MATL( I ) .EQ. ’AIRSPACE’ ) THEN 

IF  ((L(I).LE.O. .AND.I.NE.N) .OR. (HEIGHT(I) .LE.0. .AND.I.NE.N) 

2 .OR. (HEIGHT(I) .NE. 0. .AND. I .EQ.N) .OR.EMMIS(I) .LE. 0. ) THEN 

WRITE  (6, ’(IX, A, 13)’)  ’INCORRECTLY  SPECIFIED  AIR  SPACE.’, I 
STOP  ’INPUT  DATA  ERRORS,  AIR  SPACE’ 

END  IF 
ELSE 

IF  (L(I ) .LE.0 . .QR.HEIGHT(I) .LE. 0. .OR.WIDTH(I) .LE. 0. .OR. 

2 K(I ) . LE . 0 . . OR.Et“MIS(I ) .LE.0.)  THEN 

WRITE  (6, ’(IX, A, 13)’)  'INCORRECTLY  SPECIFIED  PROTECTOR I 
STOP  ’INPUT  DATA  ERRORS,  SOLID  PROTECTOR' 

END  IF 
END  IF 

40  CONTINUE 

IF  (XWALL.LT. 0. .OR. YWALL.LT. 0. ) THEN 
WRITE  (6,*)  ’INCORRECTLY  SPECIFIED  POINT  ON  WALL.’ 

STOP  'DATA  INPUT  ERRORS,  XWALL  & YWALL’ 

END  IF 

IF  (T(O).LE.T(N))  THEN 

WRITE  (6,*)  'INCORRECTLY  SPECIFIED  ENDPOINT  TEMPERATURES.’ 

STOP  'DATA  INPUT  ERRORS,  T(0)  & T(N) ’ 

END  IF 

IF  (T(O).LE.O. .OR.WIDTH(O). LE.0. .OR. HEIGHT(O) .LE.0. .OR.EmiS(O) 

2 .LE.0.)  THEN 

WRITE  (6,*)  ’INCORRECTLY  SPECIFIED  STOVE  PARAMETERS.' 

STOP  'DATA  INPUT  ERRORS,  STOVE’ 

END  IF 
END  IF 

IF  (MATL(N).EQ. ’AIRSPACE’ ) HEIGHT (N)=HEIGHT(N-1) 
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224 

435 

jc 

1 

225 

436 

1 

IF  (IVAR.EQ.O)  THEN 
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PRINT  OUT  A DESCRIPTION  OF  THE  STOVE  / WALL  SYSTEM  AS  SPECIFIED 
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WRITE  (6,11)  TITLE, N+l.XWALL, YWALL 
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IF  (IVAR.GT.0. AND.IVAR.LT. 4)  THEN 

1 
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WRITE  (6,5)  VAR ( IVAR ) , VLOWER , VUPPER , VINCR 

1 
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ELSE  IF  (IVAR.GT.4)  THEN 
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WRITE  (6,6)  VAR(IVAR) ,ISUB, VLOWER, VUPPER, VINCR 
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1 

END  IF 
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WRITE  (6,12) 
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DO  50  1=0, N 
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IF  (I.EQ.O.OR.I.EQ.N)  THEN 
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IF  (MATL(I).EQ. ’STOVE’)  THEN 
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WRITE  (6,13)  I,MATL(I),HEIGHT(I),WIDTH(I),EMlIS(I),T(I)-273. 
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ELSE  IF  ( MATL( I ).EQ. ’AIRSPACE’)  THEN 

246 

457 

1 

WRITE  (6,14)  I,MATL(I),L(I),Ef*JIS(I),T(I)-273.15 
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WRITE  (6,15)  I ,MATL(I ) , HEIGHT(I ) , WIDTH( I)  ,L(I ) , EMUS ( I ) ,K(I ) 
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IF  (MATL(I).EQ. 'STOVE')  THEN 
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WRITE  (6,13)  I ,MATL(I ) , HEIGHT (I ) , WIDTH( I)  ,EM1IS(I ) 
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ELSE  IF  ( MATL( I ).EQ. ’AIRSPACE’)  THEN 
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WRITE  (6,14)  I,MATL(I),L(I),EmiS(I) 
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WRITE  (6,15)  I ,MATL(I ) , HEIGHT(I ) , WIDTH( I) ,L(I ) , EMMIS(I) ,K(I ) 
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FORMAT  (’OSYNTAX  ERROR  ON  "FOR”  STATEMENT,  ILLEGAL  OR  NO  VARIABL| 
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2E  SPECIFIED.’,/,’  ’,A79,//) 
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FORMAT  ('OSYNTAX  ERROR  ON  ’’FOR’’  STATEMENT,  SUBSCRIPT  REQUIRED 

FO| 
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2R  VARIABLE  SPECIFIED.’,/,’  \A79,//) 

1 
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FORMAT  (’OSYNTAX  ERROR  ON  STATEMENT,  UNRECOGNIZED  KEYWORD.',/, 

1 
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2 ’ ’ , A79, //) 
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FORMAT  (’OSYNTAX  ERROR  ON  "FOR"  STATEMENT,  RANGE  AND  INCREMENT 
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2NCONSISTENT .’,/,’  \A79,//) 
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FORMAT  ( ’ 0SERIES  OF  CALCULATIONS  FOR  VARIABLE  ’ ,A, ’ : ’.3F10.3) 
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FORMAT  < ’ 0SERIES  OF  CALCULATIONS  FOR  VARIABLE  \A,’(’,I2, 
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FORMAT  (’1’ ,A, // , 
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2 ’ NUMBER  OF  NODES  IN  CALCULATION:  ’,12,/,’  POINT  ON  WALL:  (X) : 
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3 F10.3,T41, ’ ( Y ) : ' ,F10.3) 

1 

277 

488 

1 12 

FORMAT  (’01  MATERIAL  HEIGHT  WIDTH  THICK  EMISS| 
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491  | 4 (W/m  K)  CC)’,/,’  ,,7BV~,),n  I 281 

492  j 13  FORMAT  ( IX,  12 , 2X, A, T26,  IX, F6 . 2, T35, F6 .2 ,T51  ,F6 . 2, T70 , F7 . 2 ) j 282 

493  | 14  FORMAT  ( IX, 12 , 2X, A, T43,F6 .2, T51 ,F6. 2, T70 , F7 . 2)  | 283 

494  j 15  FORMAT  ( IX. 12 , 2X, A, T26, IX, 4 (F6. 2, 2X) , F8 . 3 , 2X, F8 .2)  j 284 

495  | END  | 285 
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DOUBLE  PRECISION  FUNCTION  FVAL  (IN.ICHR) 
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ARGUMENTS : 

IN: 

STRING  CONTAINING  (MAYBE)  NUMBER 

c 

ICHR: 

(INPUT)  BEGINNING  CHARACTER  POSITION 

c 

(OUTPUT)  NEXT 

CHARACTER  POSITION 

C 

C PURPOSE:  DECODE  NEXT  NUMBER  IN  STRING  AS  A DOUBLE  PRECISION  VALUE 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

CHARACTER  IN* (* ) , FORMAT*10 
IL=LEN(IN) 

IFIRST=ICHR 
DO  20  I=IFIRST, IL 

IF  ( (IN(I : I ) .GE . ’ 0’ . AND . IN( I : I) .LE. ’ 9’ ) .OR. IN(I : I ) . EQ. ’ . ' 

2 .OR. IN(I : I ) .EQ. ’+’ .OR. IN(I : I ) . EQ. ' ) THEN 
C 

C THERE  IS  A NUMBER  ON  THE  CARD,  FIND  OUT  WHAT  IT  IS 
C 

DO  10  J=I , IL 
C 

C IF  WE  FIND  THE  END  OF  THE  NUMBER,  READ  IT  FROM  THE  LINE 
C 

IF  ( ( INC J : J) .LT. ’ 0’ .OR. INC J: J) .GT . ' 9’ ) . AND. INC J: J) . NE . ’ . ’ 

2 .AND.INCJ: J) .NE. ’+’ . AND. INC J: J) . NE . ’ ) THEN 
WRITE  (FORMAT, 30)  J-I 
READ  (IN(I : J- 1 > .FORMAT)  VAL 
FVAL=VAL 
ICHR=J 
RETURN 
END  IF 

10  CONTINUE 

C 

C IF  WE  GET  TO  THE  END  OF  THE  LINE  WITHOUT  FINDING  END  OF  NUMBER, 

C JUST  READ  THE  NUMBER 
C 

WRITE  (FORMAT, 30)  IL-I+1 
READ  (INC I :IL), FORMAT)  VAL 
FVAL=VAL 
ICHR=J 
RETURN 
END  IF 

20  CONTINUE 

C 

C IF  NO  NUMBER  IS  ON  THE  CARD,  JUST  RETURN  A 0. 

C 

FVAL=0 . 

RETURN 

C 

30  FORMAT  (’ (F’,12. 2, ’ .0)’ ) 

END 
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SUBROUTINE  HEATS 


C 


c 

@@@@ 

8 

c 

@888 

8 

8 

c 

@888 

8 

@88888 

@88888 

888888 

@88888 

c 

@8888888 

8 @88 

8 

8888 

888 

c 

@880 

8 

@88888 

@88888 

@888 

@88888 

c 

8888 

8 

8 

8 @88 

8888 

888 

c 

@888 

8 

@88888 

@88888 

8888 

@88888 

c 

C ARGUMENTS : NONE 

C 

C PURPOSE:  CALCULATES  TOTAL  HEAT  AND  TOTAL  RESISTANCE  THROUGH  WALL 

C PROTECTORS 

C 

INCLUDE  ’ STOVE. CMN’ 

C 

C JUST  SUM  UP  RESISTANCES  TO  MAKE  UP  TOTAL  RESISTANCE 
C 

RTOTAL=0 . 

DO  10  1=1, N 
CALL  HEAT (I) 

1 0 RTOTAL=RTOTAL+R ( I ) 

C 

C TOTAL  HEAT  IS  JUST  DELTA  T / TOTAL  RESISTANCE 
C 

QTOTAL= ( T ( 0 ) - T ( N ) ) /RTOTAL 

RETURN 

END 
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597 
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10 

11 
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21 
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28 

29 

30 

31 

32 

33 
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35 
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37 

38 

39 

40 
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SUBROUTINE  HEAT ( I ) 


C 


c 

sees 

e 

c 

sees 

@ 

e 

c 

sees 

@ 

eeeeee 

eeeess 

c 

@@@@@@@@ 

@ eee 

@ 

sees 

c 

ssse 

e 

eeeeee 

eeeeee 

sees 

c 

sees 

§ 

§ 

e ess 

esse 

c 

esse 

@ 

seeeee 

essees 

esse 

C 


C ARGUMENTS:  I:  ELEMENT  NUMBER 

C 

C PURPOSE:  CALCULATES  HEAT  THROUGH  ELEMENT  NUMBER  I 

C 

INCLUDE  ’ STOVE. CMN’ 

DATA  SIGMA  / 5.67D-8  / 

IF  (MATL(I) .NE. ’AIRSPACE’ ) THEN 
C 

C MATERIAL  IS  A SOLID,  CONDUCTION  ONLY 
C 

R(I)=L(I)/K(I) 

ELSE  IF  (MATL(I).EQ. ’AIRSPACE' .AND. I. EQ.l)  THEN 
C 

C MATERIAL  IS  AND  AIRSPACE  NEXT  TO  STOVE,  CALCULATE  FOR  POINT  ON  WALL 
C 

XNOTS=XWALL- (XSTOVE-WSTOVE* . 5 ) 

YNOTS=YWALL - YSTOVE 
ZNOTS=ZSTOVE 

CF ( I ) =CFIGS (XNOTS , YNOTS , ZNOTS ) -CFIGS (XNOTS-WSTOVE , YNOTS , ZNOTS ) 

2 +CFIGS (XNOTS-WSTOVE , YNOTS-HSTOVE , ZNOTS ) 

3 -CFIGS (XNOTS, YNOTS-HSTOVE, ZNOTS) 

RAD (I )=1 . / (CF(I )*SIGMA*HHIS(I )*(T(I-1)**2+T(I)**2) 

2 *(T(I-1)+T(I))) 

RMTEMP=293 . 15 

RLOSS=  ( T ( I - 1 ) -T  ( I ) ) / ( ( T ( I ) **4 -RMTEMP**4 ) *SIGMA*QHIS  ( I ) 

2 *(1-CF(I) ) ) 

R(I)=1./(1. /RAD(I)-1. /RLOSS) 

ELSE  IF  (MATL(I).EQ. ’AIRSPACE’ .AND. I. NE.N)  THEN 
C 

C MATERIAL  IS  AN  AIRSPACE,  BUT  NOT  NEXT  TO  STOVE 
C 

CF(I)=CFIGP(HEIGHT(I) ,WIDTH(I ) ,L(I) ) 

H(I )=HCONV(T( 1-1) , T(I ) , HEIGHT (I ) ) 

CONV(I)=l./H(I) 

RAD(I )=1 . / (CF (I )*SIGMA*EMMIS( I )*(T( I-l)**2+T( I )**2) 

2 *(T(I-1)+T(I))) 

R(I)=l./(l./CONV(I)+l./RAD(I)) 

ELSE 

C 

C MATERIAL  IS  AN  AIRSPACE  AND  LAST  ELEMENT  ...  THE  GREAT  OUTDOORS 
C 

CF(I)=1 . 

H(I )=HCONV(T( 1-1) , T(I ) , HEIGHT (I ) ) 

CONV(I)=l./H(I) 

RAD(I)=1./(CF(I)*SIGMA"EEWIS(I)*(T(I-1)**2+T(I)**2) 

2 *(T(I-1)+T(I))) 

R(I)=l./(l./CONV(I)+l./RAD(I)) 

END  IF 

Q(I)=(T(I-1)-T(I))/R(I) 

RETURN 

END 
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644 

1 

DOUBLE  PRECISION 

FUNCTION  CFIGS  (A,B,C) 

1 

1 

645 

|c 

1 

2 

646 

|c 

eeeeeeee 

8880 

0800 

1 

3 

647 

|c 

8 8880 

8 

8 

4 

646 

ic 

0 8008 

888888  0888 

880800 

688880 

5 

649 

jc 

8 

8888 

8888 

000880 

868 

1 

6 

650 

jc 

8 8888 

8888 

8880 

8 

668088 

1 

7 

651 

|c 

8 8888 

8088 

0880 

88888080 

680 

8 

652 

ic 

08888800 

8888 

0008 

8 8860 

866888 

9 

653 

ic 

80066688 

10 

654 

ic 

1 

11 

655 

ic 

ARGUMENTS: 

A: 

WIDTH  OF  RECTANGLE 

12 

656 

ic 

B: 

HEIGHT  OF  RECTANGLE 

1 

13 

657 

|c 

C: 

DISTANCE  TO  POINT  OF  CALCULATION 

14 

658 

ic 

1 

15 

659 

ic 

PURPOSE: 

CALCULATES  RADIATION  CONFIGURATION  FACTOR  FOR  A PLANE 

16 

660 

|c 

ELEMENT 

TO  A PLANE 

PARALLEL  RECTANGLE. 

1 

17 

661 

ic 

1 

18 

662 

ic 

SOURCE:  THERMAL  RADIATION  HEAT  TRANSFER, 

SEIGEL  & HOWELL. 

19 

663 

ic 

1 

20 

664 

1 

IMPLICIT 

DOUBLE 

PRECISION 

(A-H.O-Z) 

1 

21 

665 

PI=3. 14159 

22 

666 

1 

X=A/C 

1 

23 

667 

Y=B/C 

24 

668 

CFIGS=1 . 

/(2*PI)* 

(X/ SQRT ( 1+X**2) *ATAN(Y/SQRT ( 1+X**2 ) ) 

25 

669 

1 

2 + Y/SQRT(1+Y**2)*ATAN(X/SQRT(1+Y**2) ) ) 

26 

670 

RETURN 

27 

671 

1 

END 

1 

28 

50 


672 

1 

DOUBLE  PRECISION  FUNCTION  CFIGP  (A,B,C) 

1 

673 

|C 

2 

674 

jc 

@@@@@@@@ 

mm  mm 

3 

675 

\C 

@ @@@@ 

@ 

@ 

4 

676 

jc 

@ @@@@ 

mmm  mm 

@@@@@@ 

@@@@@@  i 

5 

677 

|c 

@ 

mm  mm 

@m  e 

6 

678 

jc 

@ @@@@ 

mm  mm 

@ 

@@@  e 

7 

679 

jc 

@ @@@@ 

mm  mm 

@@@@@@@@ 

@@@@@@  i 

8 

680 

jc 

@@@@@@@@ 

mm  mm 

@ @@@@ 

@g@  ! 

9 

681 

jc 

@@@@@@@@ 

@@@ 

10 

682 

jc 

11 

683 

jc 

ARGUMENTS : 

A:  WIDTH  OF  RECTANGLE 

12 

684 

jc 

B:  HEIGHT  OF  RECTANGLE 

13 

685 

jc 

C:  DISTANCE  BETWEEN  RECTANGLES 

14 

686 

jc 

15 

687 

jc 

PURPOSE:  CALCULATES  RADIATION 

CONFIGURATION  FACTOR  FOR  TWO 

16 

688 

jc 

IDENTICAL,  PARALLEL, 

DIRECTLY  OPPOSED  RECTANGLES. 

17 

689 

jc 

18 

690 

jc 

SOURCE:  THERMAL  RADIATION  HEAT  TRANSFER, 

SEIGEL  & HOWELL. 

19 

691 

jc 

20 

692 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

21 

693 

1 

PI=3. 14159 

22 

694 

X=A/C 

23 

695 

1 

Y=B/C 

24 

696 

1 

CFIGP=2 . / (PI*X*Y)*(LOG( ( ( 1+X*X)*( 1+Y*Y) )/ ( 1+X*X+Y*Y) )**0 . 5 + 

25 

697 

2 X*SQRT(1+Y*Y)*ATAN(X/SQRT(1+Y*Y) ) + 

26 

698 

1 

3 Y*SQRT(1+X*X)*ATAN(Y/SQRT(1+X*X) ) - 

27 

699 

1 

4 X*ATAN(X) 

- Y*ATAN(Y)) 

28 

700 

1 

RETURN 

29 

701 

1 

END 

30 

51 


702 

1 

DOUBLE  PRECISION  FUNCTION  HCONV(Tl,T2,L) 

1 

1 

703 

|C 

1 

2 

70*. 

|c 

gggg  g 

1 

3 

705 

|c 

gggg  g 

4 

706 

ic 

gggg  g gggggg  gggggg  gggggg 

ggg  g 

1 

5 

707 

ic 

gggggggg  g ggg  ggg  g g ggg 

ggg  g 

1 

6 

708 

ic 

gggg  g g ggg  g g ggg 

ggg  g 

1 

7 

709 

ic 

gggg  g g ggg  ggg  g g ggg 

gggg 

8 

710 

ic 

gggg  g gggggg  gggggg  g ggg 

gg 

1 

9 

711 

ic 

10 

712 

ic 

ARGUMENTS:  Tl:  TEMPERATURE  OF  HOTTER  SURFACE 

1 

11 

713 

ic 

T2:  TEMPERATURE  OF  COOLER  SURFACE 

1 

12 

714 

ic 

L:  HEIGHT  OF  SURFACES 

13 

715 

ic 

1 

14 

716 

ic 

PURPOSE:  CALCULATES  FREE  CONVECTION  HEAT 

TRANSFER  COEFFICIENT  FOR  j 

15 

717 

ic 

A VERTICAL  SURFACE. 

1 

16 

718 

ic 

1 

17 

719 

ic 

SOURCE:  FUNDAMENTALS  OF  HEAT  TRANSFER,  INCORPERA  & DEWITT.  ! 

18 

720 

ic 

1 

19 

721 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

20 

722 

1 

DOUBLE  PRECISION  K,NU,L,NUSELT 

21 

723 

1 

DIMENSION  C(3 , 5) 

1 

22 

724 

1 

DATA  ((C(I, J),J=1,5),I=1,3)  / 

23 

725 

1 

2 -. 381021E-2,  . 132063E-3 , -.117332E-6, 

.687499E-10,  - 

. 127680E-13 | 

24 

726 

3 , - . 167333E-4 , .143076E-6,  - . 249135E-10 

.781850E-13,  - 

. 127693E-16 | 

25 

727 

4 , - . 501195E-5,  . 468550E-7 , .881329E-10, 

- . 117315E-13 , 

. 307192E-17  j 

26 

728 

1 

5 / 

1 

27 

729 

1 

TF=(Tl+T2)/2 

28 

730 

1 

DELTAT=(T1-T2) 

1 

29 

731 

1 

K=C ( 1 , 1 )+C ( 1 , 2 ) *TF+C (1,3) *TF*TF+C ( 1 , 4 ) *TF**  3+C (1,5) *TF**4 

30 

732 

ALPHA=C  ( 2 , 1 )+C ( 2 , 2 ) *TF+C ( 2 , 3 ) *TF*TF+C ( 2 , 

4)*TF**3+C(2, 5)*TF**4  i 

31 

733 

1 

NU=C(3 , 1)+C(3 ,2)*TF+C(3 , 3 )*TF*TF+C(3 , 4)*TF**3+C(3 , 5)*TF 

**4  1 

32 

734 

1 

PR=NU/ ALPHA 

33 

735 

RA=9 . 8* ( 1 . / TF ) *ABS ( DELTAT ) *L**3 / ( NU*ALPHA ) 

34 

736 

1 

NUSELT=(0. 825+0. 387*RA**(l./6. ) / ( 1. +(0. 492/PR)**(9. /16 . 

)) 

35 

737 

1 

2 **(8./27. ))**2 

36 

738 

1 

HCONV=NUSELT*K/L 

37 

739 

RETURN 

1 

38 

740 

1 

END 

1 

39 

741 

1 

CHARACTER* 5 FUNCTION  CMET (VALUE , EPS ) 

1 

1 
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|C 
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jc 
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|c 

g gggg 
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745 

ic 
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ic 

g 
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g gggg 
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7 

748 

ic 
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ic 
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g g ggg  gggggg  gggg 

9 

750 

ic 

10 

751 

ic 

ARGUMENTS: 

VALUE:  NUMBER  TO  BE  EVALUATED 

11 

752 

ic 

EPS : ACCEPTANCE  CRITERION  FOR  VALUE 

12 

753 

ic 

13 

754 

ic 

PURPOSE:  FUNCTIONS  RETURNS  A CHARACTER  INDICATION  OF  WHETHER  THE 

14 

755 

ic 

VALUE  IS  WITHIN  LIMITS.  USED  FOR  DEBUG  PRINTOUT 

15 

756 

ic 

16 

757 

IMPLICIT 

DOUBLE  PRECISION  (A-H,0-Z) 

1 

17 

758 

i 

IF  (ABS (VALUE ) . LE . EPS ) THEN 

18 

759 

CMET=’ 

(IN)  ' 

19 

760 

i 

ELSE 

20 

761 

CMET=' (OUT) ’ 

21 

762 

i 

END  IF 

22 

763 

i 

RETURN 

23 

764 

i 

END 

1 

24 
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765 

766  |C 


767  jc 

88888088 

768  jc 

gggg  § 

763  jc 

e§ eg  § 

8B8BB8 

gee 

§ 

ggggggg 

770  jc 

ggg  §g  ggg 
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gee 

§ 

g g g gg 

771  If* 

fix  |L 

§ee§ 

000000 

ggg 

g 

§ g ggg 

772  jc 
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§:§§ 

gee 

g 

g g ggg 

773  jc 

gggg 
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gggggg 

g g g gg 

"5  C A3SO*STS:  : SURFACE  BOBER  PCS  SCM4ARY  PSISTOUT 

776  C 


777 

778 
773 
76C 
781 
76 2 
783 
76* 
7 85 

786 

787 
7 86 
789 
79C 
7s: 
732 
753 
73* 
735 
~S€ 


PURPOSE : rSETTS  XT  ALL  TEE  VARIAELZS  3 THE  CMOS  SLOCK.  D5E0 

fce  debog  mncm. 

DdLMS  * STOPS.  CMS* 

'•SHE  (6,13)  I ,ISMLCI) , 1(1)  ,S(I) ,L(I) ,£CI),H(I) ,CFCD .2M-CS Cl) , 

• cgbv; : : . sad c ; . seight (i ; .kebthci ) ,xstovs , yscots , zstote . 

3 KSTC7VE,BSITCWS,»<ALL  VKALL . KJOIAL , QTCIAL . K 
RECURS 


FCEKfeC  . CCURSZ5T  VAXES  ECS  SPACE  M2.  MATERIAL  IS  A . . 


2 

3 

A 


5 

7 

8 
3 


ISKFSRAXRE  : . FLO . 3 . 1*1,  RSSISIASCE  : ; ,F10.3.  . 
lEICKSESS : .FID.  3.1*1.  3HE93&L  CCSiOXTiWiTX : * ,512.5,  !, 

E F10.3.T61,  ’ CXSFIGQRATTCE  ? ACTOR: 1 ,310.3, 

emissive:?:  ,?ic.3,:*i, ■ _vlividual  ezai:  ei:.3,  . 

COST . R3STS1WCS : ' , E1C . 3 , , ’SAD.  RESISTANCE : EM. 3.  , 


ESISEC  : 
XSTC7S: 
-STMTE: 
XnAX: 
ETCIAL : 


.fio.3,:*i,  ’tmnB:  -•  ,no.3,/, 

.?10.3,T31,'XSTOBE: ' ,E  II.  3 MSI.  MSXTE:  .E1D.3,  . 
. E10 . 3 , 131 , HSTOVE : ’ ,F1D.3,/, 

.E10.3.T31. TnALL:  MF1C.3,/, 

.F3II.3,T31,  'QT0T6L:’  .FIO.3,161,  5: 
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817 
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820 
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830 
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832 
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10 

11 

12 
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19 
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CHARACTER*60  FUNCTION  PMATL (MATL1 ,MATL2 ) 
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@@@@@666 

886 
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eeee  8 
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6888  8 

6686866 

888688 

888 

c 

@@@86666 

8 6 886 
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8888 
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8666 
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8868 

8 8 668 

8 888 
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886 

c 
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8 6 @@@ 

888886 

8686 

888 

ARGUMENTS:  MATL1:  DESCRIPTION  OF  MATERIAL  CLOSER  TO  APPLIANCE 

MATL2 : DESCRIPTION  OF  MATERIAL  FURTHER  AWAY  FROM  APPLIA 

PURPOSE:  CREATES  A CONCATENATED  DESCRIPTION  OF  TWO  MATERIALS  FOR 

PRINTOUT 


NCE 


CHARACTER* 40  MATL1.MATL2 
PMATL=MATL1 


DO  10  1=40, 1,-1 
IF  (MATL1 (I : I ) . NE . ' ’)  THEN 
PMATL (1+1: 1+3)=’  / ’ 

GO  TO  20 
END  IF 

10  CONTINUE 

PMATL=’  ’ 

20  DO  30  1=40, 1,-1 

IF  (PMATL(I : I ) . NE . ’ ’)  THEN 
PMATL (1+2 : 60)=MATL2 
RETURN 
END  IF 

30  CONTINUE 

PMATL =MATL2 

RETURN 

END 


o o 
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CCt-tON  BLOCK  FOR  PROGRAM  STOVE 


1 | IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z)  I 1 

2 |C  j 2 

3 j PARAMETER  (MAXPRO=20)  | 3 

* |C  | 4 

5 | DOUBLE  PRECISION  L,K  j 5 

6 | CHARACTER  MATL*40 ,VNAME*10 , TITLE*80  j 6 

7 j LOGICAL  DEBUG, FULPRT  | 7 

8 j COMMON  /NSTOVE/  T(0 iMAXFRO) ,R(0 :MAXFRO) ,L(0 :MAXFRO) ,K(0 :MAXPRO) , | 8 

9 | 2 H(0 :MAXPRO) , CF(0 :MAXPRO) ,EIWIS(0 :MAXPRO) ,Q(0 :MAXPRO) , | 9 

10  | 3 CONV(0 :MAXPRO) ,RAD(0  rMAXPRO) ,HEIGHT(0:MAXPRO) ,WIDTH(0:MAXPRO) , | 10 

11  j 4 XSTOVE , Y STOVE , ZSTOVE , WSTOVE , HSTOVE , XWALL , YWALL , R TOTAL , QTOTAL , N , j 11 

12  | 5 IVAR, ISUB,VLCWER,VUPPER,VINCR,DEBUG, FULPRT  | 12 

13  j COMMON  /CSTOVE/  MATL ( 0 : MAXPRO ) , VNAME , TITLE  | 13 

14  |C  | 14 


Appendix  fi: 


Data  Input  for  STOVE 
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The  data  input  for  stove  takes  the  form  of  six  different  key  words  with 
arguments  to  specify  values  which  depend  upon  the  key  word.  In  most  cases, 
the  order  of  the  key  words  is  unimportant,  except  as  noted  below.  A descrip- 
tion of  each  of  the  input  key  words  and  values  which  go  on  the  same  line  are 
presented  below: 


STOVE 

AIRSPACE 

FOR 

XWALL 

YWALL 

PROTECTOR 


<height>  <width>  <emissivity>  <temperature> 

<thickness>  <emissivity> 

<variable>  — <lower>  <upper>  <increment> 

<x  position> 

<y  position> 

<thickness>  <height>  <width>  <k>  <emissivity>  <temperature> 


BOLDFACE  type  are  required  key  words.  Words  in  <brackets>  specify  numeric 
inputs  as  follows : 

<emissivity>  specifies  the  emissivity  of  the  cooler  surface  of  the 

material.  If  specified  for  an  airspace,  it  is  the  emissivity 
of  the  surface  adjacent  to  the  airspace  at  the  farther 
distance  from  the  stove. 


<height>  specifies  the  height  of  the  stove  or  protector  in  meters. 

<increment>  specifies  the  amount  to  increment  the  variable  <variable>  in 

the  FOR  statement  for  each  calculation  to  be  performed.  The 
first  calculation  is  done  with  <variable>  equal  to  the  value 
<lower>;  the  second  calculation  is  done  with  <variable>  equal 
to  the  value  <lower>  + <increment>  and  so  forth  until  the 


value  of  <variable>  is  greater  than  or  equal  to  the  value  of 
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<k> 


<lower> 


<temperature> 


<thickness> 


<variable> 


<upper>.  The  units  for  the  number  are  the  same  as  those  for 
the  variable  <variable>. 

specifies  the  thermal  conductivity  of  the  solid  protector  in 
W/m-K. 


specifies  the  beginning  value  of  the  variable  <variable>  in 
the  FOR  statement  for  each  calculation  to  be  performed.  The 
first  calculation  is  done  with  <variable>  equal  to  the  value 
<lower>;  the  second  calculation  is  done  with  <variable>  equal 
to  the  value  <lower>  + <increment>  and  so  forth  until  the 
value  of  <variable>  is  greater  than  or  equal  to  the  value  of 
<upper>.  The  units  for  the  number  are  the  same  as  those  for 
the  variable  <variable>. 

specifies  the  temperature  of  the  stove  surface,  protector,  or 
airspace  in  K.  Temperatures  are  only  specified  for  the  stove 
surface  (surface  number  0)  and  for  the  outermost  surface  or 
airspace  (surface  number  N) . 

specifies  the  thickness  of  the  material  (for  a protector)  or 
the  distance  between  surfaces  (for  an  airspace) . 

specifies  the  variable  to  be  incremented  in  each  calculation 
to  be  done.  The  first  calculation  is  done  with  <variable> 
equal  to  the  value  Clowe r>;  the  second  calculation  is  done 


<width> 


<x  position> 


<y  position> 
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with  <variable>  equal  to  the  value  <lower>  + <increment>  and 
so  forth  until  the  value  of  <variable>  is  greater  than  or 
equal  to  the  value  of  <upper>.  Legal  variables  which  may  be 
used  are:  T(0)  --  the  stove  temperature,  width(i),  k(i) , 

xwall , ywall , emissivity(i) , height(i) , l(i). 

specifies  the  width  of  the  stove  or  protector  in  meters. 

specifies  the  x position  of  the  point  on  the  wall  at  which  the 
calculation  is  to  done  in  meters. 

specifies  the  y position  of  the  point  on  the  wall  at  which  the 


calculation  is  to  done  in  meters . 
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